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. Abstract 

-^ Low complexity joint estimation of synchronization impairments and channel in a single-user MIMO- 

OFDM system is presented in this letter. Based on a system model that takes into account the effects of 
synchronization impairments such as carrier frequency offset, sampling frequency offset, and symbol tim- 
ing error, and channel, a Maximum Likelihood (ML) algorithm for the joint estimation is proposed. To 
reduce the complexity of ML grid search, the number of received signal samples used for estimation need 
to be reduced. The conventional channel estimation methods using Least-Squares (LS) fail for the reduced 
sample under-determined system, which results in poor performance of the joint estimator. The proposed 

ML algorithm uses Compressed Sensing (CS) based channel estimation method in a sparse fading scenario, 
"xT ' 

where the received samples used for estimation are less than that required for an LS based estimation. The 

"xl" . 

^■i. , performance of the estimation method is studied through numerical simulations, and it is observed that CS 

based joint estimator performs better than LS based joint estimator. 
O 
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1. Introduction 

Multiple Input Multiple Output-Orthogonal Frequency Division Multiplexing (MIMO-OFDM) system, 

a ■ 

the preferred solution for the next generation wireless technologies, is very sensitive to synchronization 
impairments such as Carrier Frequency Offset (CFO), Sampling Frequency Offset (SFO) and Symbol Timing 



Error (STE) 



H- 



04|]. In this letter, we propose a low complexity Maximum Likelihood (ML) algorithm for 



the joint estimation of synchronization impairments and channel using Compressed Sensing (CS) technique, 
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in a sparse fading scenario, where the received samples used for estimation are less than that required for a 
Least Squares (LS) based estimation. 

2. System Model 

Consider a MIMO-OFDM system with Nj transmit antennas and Nr receive antennas using Quaternary 

Phase Shift Keying (QPSK) modulation and N subcarriers per antenna. Let T be the sampling time at the 
transmitter and f c be the carrier frequency. We define the normalized CFO as e = Af c NT, the normalized 
SFO as r] — AT/T, and the normalized STE as 6, where Af c is the net CFO in the received signal and AT is 
the difference between the sampling time at the receiver and the transmitter [4]. Let X be the block diagonal 
matrix with each diagonal matrix having the signal vector transmitted from each transmit antenna. Also, let 
h be the column vector representing the MIMO channel with L m as the maximum length of channel between 
any transmit and receive antenna pair. The signal vector at the receiver side is derived in |4J as, 

r = Ai(e,?7,0)h + w (1) 

where Ai(e, ?/, 9) = I Nr ® (D(e, rj)¥i(ri)G(e)X(I NT e F 2 )) 

D(e, if) = diag\\, exp(j2^e(l + rj)/N), ..., exp(j2jis(l + if)(N - l)/N)] 

G(0) = diag[l,expi-j2ne/N), ..., exp(-j2n(N - 1)8/N)], 

rp / m exp(;2^(n(l + i?))/AQ 

[Fi(T])]n,k = , [F 2 ]i,/ = exp(-j27Tlk/N), 

with n, k - 0, 1, . . . , N — 1, and / = 0, 1, . . . , L m — 1. w is the additive circular Gaussian noise vector with 
mean zero and variance cr^,. Let ma x denote the maximum STE. Then the system model in (Q]) can be 
re-written as, 



r = A2(e, r])hg + w 
where A 2 = I Nr O (D(e,^)F 1 (^)X(I A , r ®F 2 e,„„)), 
and [F 2 g m Jkj = exp(-j2nlk/N), 



(2) 



with k - 0,1,..., N -1,1 = 0,1,..., (L m + max - 1), and h fl being the STE embedded MIMO channel as 

a. 



given in | 



Notations: Upper case bold letters denote matrices and lower case bold letters denote column vectors. A denotes the estimate of A. 
\m denotes an M X M identity matrix, identity matrix. A H and A^ denote conjugate transpose, and pseudo-inverse of A, respectively. 
[A]„, „ denotes the (m, n) th element of A. ® represents Kronecker product. diag[x] represents a diagonal matrix having the elements 
of x as diagonal elements. Tr(A) represents trace of A. Calligraphic letter T denotes set and T c denotes set complement. At-(A(7-)) 
denotes the column (row) sub-matrix of A formed by the columns (rows) of A listed in the set T. 



3. ML Algorithm for Joint Estimation 

The ML cost function [4] of the parameters e, 6, rj, and h, obtained from (Q} is, 

arg min J(e, 77, 6, h|r) = arg min (r - A 1 h) // (r - A^). (3) 

The multi-dimensional minimization in (f3]l gives the estimate of the parameters e, 8, 77, and h. Given the 
estimate of channel, h and he, and using the system models in (Q]) and (0, the optimization problem in (01 
reduces to a two-dimensional and one-dimensional minimization problem respectively as, 

[e,ff\ = argmin(r- A 2 h e ) H (r - A 2 h e ) = argminJi(e,?7|r,n e ), (4) 

e,V e,i7, 

[6] = argmin(r - Aih) ff (r - Aih) = arg min J2(0|r, e, fj, h). (5) 

8 8 

For the above ML algorithm to have a unique solution with the LS estimate of the channel, the number 
of received signal samples used for estimation must at least be equal to the number of unknown channel 
coefficients, i.e., MNr > L„NtNr. To have a low complexity joint estimation at the receiver we need 
to reduce the received samples used for estimation, where the ML algorithm using LS channel estimation 
(MLLS) fails. Hence we propose an ML algorithm using CS technique which performs better than MLLS 
for an under-determined MIMO-OFDM system in sparse fading channel. 



3.1. CS based channel estimation 



CS is a novel technique where a 
fewer samples than usually required 



arameter that is sparse in a transform domain can be estimated with 



51 



16] 



]. The application of CS is to recover the ^-sparse channel (A 
channel is said to be /f-sparse if it contains at most K non-zero coefficients) from MNr received signal 
samples, where MNr < L m N T N R . Using (Q} and ©, 

r« = T(r) = Ai„(e, 77, ff)h + w„ = A 2u (e, rj)h e + w„, (6) 

where T(v) is the operator which randomly selects M samples from each receive antenna given in r. Also, 
Ai u =Ai (r) and A2„=A2(7-) where T contains the indices of the MNr samples selected from r. In CS frame- 
work, r„ is called the observation vector and A, which represents either Ai„ or A2 M , is called the measurement 
matrix. 

In this letter, we use Subspace Pursuit (SP) algorithm [7] which is a popular greedy algorithm used in 
CS. In each iteration, SP identifies a Zf-dimensional space that reduces the reconstruction error of the sparse 
channel h. The steps involved are given in Algorithm 1. It has been shown theoretically that SP algorithm 
converges in finite number of steps [7]. 



Algorithm 1 Sparse channel estimation using SP Algorithm 



Inputs: A, r, and K 
I: A = AC; 

2: Initialization: k = 0, To = 0, eo 
3: repeat 



-*- Normalize columns of A using diagonal matrix C: 



k = k+ 1; 

T = Tifc_iU (indices of ^-highest magnitude components of A H et_ 

Vf = Air, v-f r = ; 



v r - a. t i, v r 

7j = indices of A"-highest magnitude components of v 

8: e k = r - A n A^r; 

9: untll(||e,|| 2 >||en|| 2 ) 
10: n = T k .i; 
11: h n =A^r,h r c=0; 
Output: T k , h (SP) = Ch 



Algorithm 2 MLSP 



Inputs: r a , T, [0 m j n , #max, #gridL[£min, fmax, egridM^min, J/max, ?7grid] 

1: for j = 6 min : egrid : <W do 



2: 


for k = r] mm : ri %M : rj max do 


3: 


Construct A 2 „(j, k); 


4: 


Obtain n< SP) ; 


5: 


Evaluate Ji Ufc|r B ,fif p) ); 


6: 


end for 


7: 


end for 



8: [cmlsp^mlsp] = argminj, Ufclr^nP ); 

j.k \ 1* I 

9: for i = 6> min : 9 glili : 6> max do 
10: Construct A hl (i, e M LSP, >7mlsp); 
11: Obtain h (! 



u(SP). 



12: Evaluate J 2 (j|r„, e ML sp, J?mlsp, fi- 



h (SP) V 



13: end for 

14: [6mlsp] = argmin 



dn J 2 (j'|r„, cmlsp, Omlsp, nf p) ); 



15: K 



j,(SP, 
< J MLSP 



Output: [9 MLSP , 6mlsp. ^mlsp, h MLS p] 



• using (|6) 

• using Algorithm 1 

• using 10 



• using (|6) 

• using Algorithm 1 

• using l[5} 



3.2. ML algorithm using SP channel estimation (MLSP) 



To obtain MLSP, the estimate of h using SP, denoted as h g ' and h (SP) , obtained from Algorithm 1 are 
used to rewrite the cost function in (@]l and (O as, Ji (e, j]\r u , h g J and J2 (f9|r„, e, 77, h (SP) j, respectively. The 
steps involved in MLSP are given in Algorithm 2. 

Remarks.The computational complexity of LS based estimation in MLLS is approximately 0((L m NTNR) y ), 
whereas that of SP based estimation in MLSP is approximately 0(MN^NTL m K) [7] which is lesser. 

4. Simulation Results and Discussions 

We considered a 2 x 2 MIMO-OFDM system having N = 128 subcarriers for each transmitter with 20 
MHz signal bandwidth. The channel coefficients are modeled as circular complex-valued Gaussian random 
variable having unit variance, and uniform power delay profile with L,„=26 and sparsity level, K=5. Also, 
the transmitted symbols belong to QPSK constellation with unit magnitude. We considered the training 
blocks having a Cyclic Prefix (CP) of length 32. The condition (L m + 6 max ) less than length of CP [4] results 
in #max=5 and \6\ < 5. The range of normalized CFO used for grid search is |e| < 0.4 with a resolution 
of 10~ 2 and that of normalized SFO is \r]\ < 5 X 10~ 3 with a resolution of 10~ 4 . The actual values of the 
impairments, e, 77, and 9 used in the simulations are 0.102, 101 ppm, and 2, respectively. 



Synchronization Impairments: e=0.1 02, rj=1 01 ppm, 0=2 




Figure 1 : CRLB and MSE for the estimation of CFO. 



The Mean Square Error (MSE) values of the estimated parameters, using MLLS and MLSP, are calcu- 
lated and are plotted in log-scale against SNR(dB), together with Cramer-Rao Lower Bound (CRLB) of the 
parameters [4], in Fig([T).- Fig®. MLLS is simulated using MLSP algorithm given in Algorithm 2 by re- 
placing the SP estimate of channel obtained in step 4 and step 5 using LS estimate of the channel. It is found 
from Fig©.- Fig©, that the MSE plots of MLLS for the estimation of CFO, SFO, and channel for M=45 



Synchronization Impairments: £=0.1 02, ri=1 01 ppm, 6=2 
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Figure 2: CRLB and MSE for the estimation of SFO. 
Synchronization Impairments: £=0.1 02, n=1 01 ppm, 6=2 
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Figure 3: Tr(CRLB(h)) and MSE for the estimation of channel. 



fail, due to the poor performance of LS based estimation in under-determined system. Also, the MSE plots 
of MLSP for the estimation of CFO, SFO, and channel follow CRLB(e), CRLB(7/), and Tr(CRLB(h)) |4fl, 
respectively, but with a performance degradation of around 12 dB, 13 dB, and 15 dB SNRs, respectively, at 
high SNR. The Probability of Timing Failure [4] for the estimation of 6, defined as P t f (p) — Pr \\6 - 0\ > p , 
is calculated for p-2 and is plotted in Fig©, for MLLS and MLSP, respectively. As in the cases of CFO, 
SFO, and channel, MLSP performs better than MLLS for the estimation of STE also. It is observed from 
the figures that, to have a comparable performance with MLSP using 90 samples (M=45), MLLS requires 

at least 150 samples (M=75), which shows the difference in computational complexity. 
5. Conclusion 

In this letter, we presented a low complexity ML joint estimation algorithm for single-user MIMO- 
OFDM system, where the received samples used for estimation are less than that required for an LS based 
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Figure 4: Probability of Timing Failure as a function of SNR(dB). 



ML estimation, MLLS. An ML algorithm for the joint estimation of synchronization impairments and chan- 
nel using CS based technique, MLSP, is proposed. It is found from the simulations that MLSP performs 
better than MLLS for the joint estimation of CFO, SFO, STE, and channel. 

References 

[ 1 ] Morelli M, Kuo C-CJ, Pun M-O. Synchronization Techniques for Orthogonal Frequency Division Mul- 
tiple Access (OFDMA): A Tutorial Review. Proc. IEEE 2007;95(7): 1394-427. 

[2] Nguyen-Le H, Le-Ngoc T, Ko C.C. Joint Channel Estimation and Synchronization for MIMO OFDM in 
the Presence of Carrier and Sampling Frequency Offsets. IEEE Trans. Veh. Technol. 2009;58(6):3075 

[3] Jose R, K.V.S. Hari. Joint Estimation of Synchronization Impairments in MIMO-OFDM System. Pro- 
ceedings of National Conference on Communication, India, 2012.p.l-5. 

[4] Jose R, K.V.S. Hari. Maximum Likelihood Algorithms for Joint Estimation of Synchronization Impair- 
ments and Channel in MIMO-OFDM Svstem. [arXiv:1210.5314 v2[cs.ITl. 

[5] Donoho. Compressed Sensing. IEEE Trans. Inf. Theory 2006;52(4): 1289 -306. 

[6] Candes EJ, Wakin MB. An Introduction to Compressive Sampling. IEEE Signal Processing Mag. 
2008;25(2):21-30. 

[7] Wei D, Milenkovic O. Subspace Pursuit for Compressive Sensing Signal Reconstruction. IEEE Trans. 
Inf. Theory 2009;59(5):2230-49. 



